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SECTION  1 


INTRODUCTION 

This  report  presents  the  development  of  a  new  technique  for 
enhancing  the  fidelity  of  flight  simulators.  The  name  given  to  the  tech¬ 
nique  and  to  the  system  which  applies  the  technique  is  BIOCONAID,  an 
acronym  for  Bionic  Control  of  Acceleration  Induced  Dimming. 

When  an  aircraft  is  subjected  to  a  change  in  its  velocity  vector  an 
acceleration  is  experienced  by  the  occupants.  The  common  situation  is 
that  of  a  pilot  changing  the  direction  of  motion  of  the  aircraft  through  a 
steep  turn  or  pull-up.  Then,  if  the  maneuver  is  coordinated,  the  pilot 
experiences  an  acceleration  tending  to  pull  him  tighter  against  his  seat 
as  if  he  weighed  more.  This  acceleration  also  affects  the  blood  pres¬ 
sure  distribution  in  his  body.  The  blood  pressure  supplying  the  eyes  and 
the  brain  is  reduced  and  the  pilot  may  experience  gray-out  or  black¬ 
out  if  the  turn  is  continued  or  the  pull-up  becomes  steeper. 

In  a  combat  situation,  the  pilot  needs  to  continue  the  turn  or  pull-up 
and  retain  his  vision.  He  may  perform  straining  maneuvers  (described  in 
Appendix  C)  to  increase  the  eye-level  blood  pressure  and  fight  the  black¬ 
out  tendency.  Such  straining  represents  expenditure  of  energy  which 
leads  to  fatigue.  The  "G-suit"  is  a  protective  garment  which  exerts  pres¬ 
sure  on  the  lower  extremities,  thereby  restricting  the  pooling  of  blood  to 
the  lower  extremities  during  high-acceleration  maneuvers.  Pilots  of 
high  performance  aircraft  wear  this  garment  to  increase  their  tolerance 
to  acceleration. 

Modern  flight  simulators,  which  are  widely  used  for  combat  crew 
training,  provide  visual  scene  dimming  as  the  simulated  G-load  increases. 
Also  provided  are  G-suits  which  realistically  react  to  the  G-load.  Missing, 
however,  is  the  pilot's  need  to  perform  the  beneficial  straining  maneuver 
with  its  attendant  energy  expenditure.  Also  missing  is  the  characteristic 
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V 


loss  of  vision  progressing  from  the  peripheral  field  inward.  The 
BIOCONAID  System  provides  the  means  to  simulate  both  of  these  missing 
elements . 

The  BIOCONAID  System  consists  of  four  integrated  software  mode.s 
plus  hardware  which  senses  and  processes  pilot  straining  information  for 
input  to  the  computer.  The  computer  software  modeis  are: 

•  Cardiovascular  Model  (predicts  effect  of  G-level  on  eye-level 
blood  pressure) 

•  G-suit  Model  (predicts  beneficial  effect  of  G-suit  at  the  input 
G-level) 

•  Straining  Model  (predicts  the  beneficial  effect  of  straining* 

•  Visual  Field  Model  (predicts  the  visual  scene  brightness  and 
the  limits  of  the  visual  field  under  acceleration) 

The  complete  system  also  includes  an  electromyographic  (EMC)  subsystem 
(signal  processor). 

Section  2  of  this  report  includes  a  discussion  of  the  development  of 
each  of  the  models  mentioned  above  plus  the  EMG  Subsystem  (EMG  strain¬ 
ing  signal).  Also  treated  is  the  integration  of  the  entire  system.  Sections 
3  and  4  cover  development  and  documentation  of  the  algorithms  and  strain¬ 
ing  experiments. 

Section  5  gives  conclusions  and  recommendations.  Briefly,  the 
BIOCONAID  System  performs  as  planned.  A  usable,  smoothed  EMG  signal 
is  produced  by  pilot  straining.  Either  the  abdomen  (external  oblique)  or 
the  back  (latissimus  dorsil  muscles  are  suitable  as  electrode  sites. 

Neither  extraneous  motions  nor  electronic  noise  seems  to  affect  the  system 
adversely.  AU  software  algorithms  perform  satisfactorily.  Full  implemen¬ 
tation  of  the  Visual  Field  Model  was  not  accomplished  during  the  contract 
period  due  to  equipment  limitations. 

Because  possible  detrimental  cardiovascular  effects  of  sustained 
straining  in  a  1-G  environment  may  exist,  it  is  recommended  that  some 
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medical  research  be  conducted  to  identify  system  limitations  or  restrictions. 

It  is  a  is  o  recommended  that  the  BIOCON  A  ID  System,  with  possible  limitations 
or  restrictions,  be  evaluated  in  an  Air  Force  simulator  with  experienced 
pilots. 

Appendix  A  gives  a  test  plan  for  system  evaluation  which  might  bp 
accomplished  in  the  F-4  simulator  at  Luke  Air  Force  Base.  Appendix  B 
gives  a  complete  listing  of  the  computer  programs  developed  during  this 
effort.  Appendix  C  gives  a  discussion  of  the  M-l/L-1  straining  maneuvers 
and  Appendix  D  describes  the  special  Tustin  solution  method  used  for  the 
cardiovascular  model. 


SECTION  2 


THEOR  V 


The  human  organism  is  a  highly  complex  system  of  interconnecting  and 
interdependent  systems.  Its  responses,  in  general,  must  be  considered  non¬ 
linear  and  extremely  difficult  to  model.  However,  bv  restricting  the  number 
and  levels  of  variables  in  a  specific  situation,  it  is  often  possible  to  model 
complex  systems  succes sfully. 

The  purpose  here  is  real-time  modelling  of  the  human  response  to  G- 
loading.  To  simplify  the  problem,  it  is  assumed  that  only  z-directed  (normal) 
acceleration  will  be  considered  and  that  the  onlv  variables  influencing  accelera¬ 
tion  are  G-level,  G-suit  protection,  and  pilot  straining.  The  complex  human 
system  can  then  be  partitioned  into  four  tractable  models:  the  cardiovascular 
model,  the  protective  garment  model,  the  straining  model  and  the  visual  scene 
model.  The  relationship  between  these  models  is  shown  in  Figure  1.  The 
result  of  each  of  these  blocks  is  a  contribution  to  the  pseudo  blood  pressure 
variable  called  "protection  variable"  (PV).  In  all  cases,  the  contribution  to 
this  blood  pressure  variable  is  either  derived  directly  from  the  acceleration 
profile  or  is  a  secondary  result  of  that  profile  and  the  operation  of  life  support 
equipment.  Note  that  a  set-point  value  (one-G  blood  pressure  at  eye  level)  is 
summed  with  the  PV  values  (see  Figure  1).  A  normally  distributed  random 
variation,  labeled  n"  is  also  included  to  account  for  variations  with  pilots  on 
different  days. 

The  cardiovascular  response  model  is  represented  bv  a  linear  transfer 
function  which  corresponds  to  the  acceleration  response  of  the  pilot's  blood 
pressure.  One  prime  governing  factor  in  pilot  response  to  acceleration  is  the 
onset  of  grayout  and  blackout.  These  visual  problems  are  directly  related  to 
the  blood  pressure  at  eye  level.  This  model  output  provides  a  dynamically 
correct  signal  which  is  equivalent  to  the  nominal  loss  in  eve  level  blood 
pressure  for  an  unprotected  human  undergoing  any  acceleration  profile,  G(t). 
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PROJECTION 


Figure  1.  System  Design  Plan 


The  straining  simulation  model  accounts  for  the  G  tolerance  enhance¬ 
ment  which  is  afforded  by  a  properly  executed  M-l  or  L-l  straining  maneu¬ 
ver.  The  purpose  of  these  straining  maneuvers  is  to  increase  the  eye  level 
blood  pressure  in  order  to  maintain  vision.  Proper  performance  of  the 
maneuvers  requires  that  the  abdominal  and  upper  torso  muscles  be  tensed 
isometrically  and  that  expirations  should  be  made  against  a  closed  glottis. 

The  result  is  an  increased  intrathoracic  pressure  and  increased  blood  pres¬ 
sure  at  the  eye.  Proper  application  of  the  straining  maneuvers  also  results 
in  the  appearance  of  myoelectric  signals  on  the  skin  surface.  These  low 
voltage  electric  signals  are  directly  associated  with  muscle  contraction 
(Bigland  and  Iippold,  1954;  Brody,  Scott,  and  Balasubiramanian,  1974). 

These  biologically  derived  signals  are  then  processed  by  the  model  to  gener¬ 
ate  a  straining  contribution  to  the  protection  variable,  PV  ,  which  represents 

s 

the  increased  blood  pressure  due  to  the  M-l  or  L-l  maneuver. 

The  pressurized  G  suit  is  a  protective  garment  used  to  increase  the 

individual's  tolerance  to  +G  .  The  suit  uses  pressurized  bladders  to  press 

z 

against  the  legs  and  lower  abdomen.  The  external  pressure  inhibits  pooling 
of  the  blood  volume  in  the  lower  extremities,  thus  insuring  a  better  supply 
to  the  heart  during  acceleration.  The  suit  must  be  inflated  by  the  G  valve 
to  a  predetermined  level  for  the  suit  to  be  effective.  The  simulation  model 
accounts  for  the  required  pressure  level  and  uses  the  actual  suit  to  provide 
the  necessary  dynamics.  The  suit  pressure  is  compared  with  the  normal 
flight  pressure  schedule  and  a  protection  variable  contribution  is  generated 
by  the  model. 

The  dynamic  visual  field  model  was  developed  to  be  readily  implemented 
in  a  simulation  system.  The  visual  field  model  uses  the  eye  level  blood  pres¬ 
sure  predicted  from  the  simulated  G  profile  and  the  other  protection  variable 
contributions  as  inputs  to  produce  dynamically  responsive  signals  which  pre¬ 
dict  the  expected  visual  field  of  a  pilot  undergoing  that  G  profile. 

The  system  presented  on  the  following  pages  represents  the  results  of 
partitioning  this  complex  physiological  system  into  models.  Each 


model  is  explored  and  developed  in  detail  based  on  a  common  protection 
variable,  F\J  which  is  related  physiologically  to  system  blood  pressure. 

The  blood  pressure  models  are  then  finally  combined  in  a  systematic 
paradigm  which  provides  the  driving  values  for  the  visual  field  response 
model. 

2.  1  THE  CARDIOVASCULAR  RESPONSE  MODEL 

Human  tolerance  to  long  term  +G  acceleration  is  normally  measured 
in  terms  of  visual  loss  (blackout)  and  unconsciousness.  Both  of  these 
tolerance  end  points  are  related  to  the  ability  of  the  cardiovascular  system 
to  deliver  oxygenated  blood  at  adequate  pressure  to  the  retinal  and  cerebral 
regions.  The  acceleration  causes  a  changing  blood  pressure  profile  in 
the  human  such  that  the  effective  pressure  vertically  above  heart  level  is 
decreased  and  the  pressure  below  heart  level  is  increased.  There  is 
therefore  a  lower  perfusion  pressure  at  eye  level.  The  distribution  of  the 
blood  in  the  body  also  changes  as  the  acceleration  pools  blood  in  the  lower 
parts  of  the  body  and  lungs.  There  is  therefore  less  available  blood  to 
circulate  and  a  lower  oxygen  content  because  the  lungs  do  not  operate  as 
efficiently. 

There  are  two  cardiovascular  systems  which  are  dynamically  involved 
in  the  process  of  blood  pressure  maintenance  while  the  human  is  under¬ 
going  +  G  acceleration.  The  hydrostatic  system  is  responsible  for  the 
z 

reduced  retinal  perfusion  pressure  at  the  eye  and  eventual  loss  of  pressure 
at  the  cerebral  level.  The  orthostatic  system  is  related  to  blood  pooling  in 
the  lower  body  with  a  concomitant  reduction  of  venous  return  to  the  heart. 

The  cardiovascular  system  has  feedback  systems  which  are  affected 
by  blood  volume  and  pressure.  The  feedback  systems  attempt  to  regulate 
the  pressure  and  flow  characteristics  of  the  cardiovascular  system.  One 
of  the  primary  pressure  sensors  for  this  system  is  located  in  the  carotid 
artery.  Any  change  in  pressure  at  the  sensors  results  in  a  regulatory 
operation  which  attempts  to  return  the  pressure  to  a  preset  value.  For 


7 


example,  a  pressure  drop  at  the  carotid  sinus  initiates  a  neurologic  reflex 

action  which  changes  the  heart  rate  and  output  so  as  to  increase  the  overall 

pressure  level.  This  section  develops  the  rationale  for  a  cardiovascular 

transfer  function  model  which  describes  the  G  acceleration  response  of  the 

z 

blood  pressure  delivery  system. 

A  radically  new  model  for  the  cardiovascular  system  was  not  the 
intent  of  this  system.  Rather  a  reasonably  accurate  transfer  function 
based  on  physiological  and  physical  principles  as  well  as  experimental 
observation  was  sought.  The  major  requirements  were  that  the  transfer 
function  model  retain  the  major  response  characteristics  of  the  cardio¬ 
vascular  system,  be  implementable  in  a  real  time  digital  simulation  system, 
and  retain  factors  which  are  identifiable  with  measurable  physiological 
factors. 

A  number  of  models  of  the  cardiovascular  system  have  been  developed 
(Gillingham,  Freeman,  and  McNee,  1972;  Knapp,  Randall,  Evans  and 
Marquis,  1978;  Koushanpour  and  Spickler,  1975;  and  Miller  and  Green, 

1973).  The  model  developed  by  Gillingham  (Gillingham,  et  al. ,  1972)  is 
used  in  the  present  effort.  This  linear  model  uses  a  transfer  function  of 
the  form: 


P(s) 


Kx(l  +  al8) 

1  +  bjS  +  b2s* 


pVc(s) 

G(s) 


(1) 


Values  selected  for  a^,  b^  and  b^  are  compromise  values  from  the 
literature  and  Gillingham's  response  curves.  The  break  frequency  for 
the  lead  compensator  factor,  (1  +  a^s),  was  selected  at  30  mHz  (millihertz). 
The  denominator  factor  (1  +  b^s  +  b^s^)  was  selected  to  give  complex  roots 
with  a  natural  frequency  of  70  mHz  and  a  damping  ratio  of  0.7.  These 
characteristics  result  in  time  responses  and  Bode  plots  which  agree 
reasonably  well  with  experimental  measurements.  To  obtain  these 
characteristics  we  find  that  a j  =  5.  31,  bj  =  3.  23  and  b£  =  5.  17. 


In  order  to  evaluate  the  system  gain  constant,  K^,  Rogers,  1978, 
treated  the  artery  from  the  aortic  arch  to  the  eye  as  a  column  of  blood 
under  the  hydrostatic  action  of  the  G-loading.  This  column  is  offset  from 
the  Z-axis  by  an  angular  displacement  of  8  degrees  due  to  the  angular 
inclination  of  the  pilot's  seat.  The  system  gain  determined  in  this  way 
=  -21.4,  compares  very  closely  with  the  gain  derived  from  experi¬ 
mental  blackout  results  (Burton,  Leverett  and  Michaelson,  1974).  When 
the  seat  inclination  is  included,  becomes  -21.4  cos  8.  The  final 
transfer  function  is  shown  in  Figure  2.  (Note  that  G(t)  here  is  dynamic 

G-loading  =  G  -1.)  The  solution  technique  used  for  generating  PV 
z  c 

from  G(t)  is  described  in  Appendix  D. 


G(s) 

-21.4  cos  0(1  +  5.  31s) 

PVc 

- > 

1  +  3.  23s  +  5. 17s2 

- > 

Figure  2.  Block  Diagram  of  Cardiovascular  System  Response 

2.2  THE  PROTECTIVE  GARMENT  MODEL  (G-SUIT  MODEL) 

The  history  of  G-  suits  dates  from  World  War  II  and  has  been  well 
reviewed  in  recent  monographs  (Burton,  et  al.,  1 974;  Fraser,  1966; 

Gaur  and  Zuidema,  1961;  and  Howard,  1965).  The  wrap-around  CSU-3/P 
cutaway-type  of  anti-G  suit,  presently  used  by  the  U.S.  Air  Force, 
improves  blackout  tolerance  in  the  +G^  direction  by  about  2G  above 
resting  tolerance  (Leverett,  Whitney  and  Zuidema,  1961).  However,  the 
mechanism  by  which  this  beneficial  effect  takes  place  is  not  entirely 
understood.  Investigations  have  shown  that  the  protection  afforded  by 
the  suit  is  not  associated  with  any  increase  in  cardiac  output  (Lindberg 
and  Wood,  1963)  nor  with  any  increase  in  eye-level  blood  pressure 
(McCally,  1970).  The  G-suit  attempts  to  prevent  pooling  of  the  blood  in 
the  lower  extremities  of  the  body  during  plus-G  maneuvers  (Leverett, 
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et  al. ,  1961).  Thus,  it  is  not  possible  to  model  the  physiological  system 
dynamics  pertaining  to  the  G-suit  until  more  knowledge  is  gained  about  the 
system. 

The  approach  taken  in  the  present  research  effort  is  to  use  an  actual 
valve  plus  a  G-suit  on  the  pilot  to  provide  the  suit  dynamics  and  to  model 
only  the  known  beneficial  effect--  the  2G  improvement  in  blackout 
tolerance  with  a  properly  operating  suit.  The  valve  will  be  activated  by 
the  computer-generat- d  G  signals. 

The  normal  suit  pressure  versus  G-loading  schedule  is  shown  in 
Figure  3.  If  the  actual  suit  pressure,  at  a  given  G-level,  falls  in  the 
cross-hatched  area  labeled  "Acceptable  Region",  then  partial  protection 
is  being  provided.  Full  2G  protection  is  afforded  when  the  suit  pressure 
reaches  the  right-hand  edge  of  the  "Acceptable  Region".  At  higher  suit 
pressures  no  additional  protection  is  provided  and  at  suit  pressures  to  the 
left  of  the  "Acceptable  Region"  no  protection  is  provided. 

A  linear  variation  of  protection  within  the  acceptable  region  is 
assumed  which  has  a  width  of  approximately  80  mm  Hg.  The  simulation 
procedure  is  to  obtain  the  difference  between  the  actual  suit  pressure  and 
the  standard  curve  pressure  (right  hand  edge  of  acceptable  region, 

Figure  3).  This  difference  is  labeled  4P  .  Then  the  protective  garment 
protection  value  contribution  is  generated  as  follows: 

if  A  P_  <  -80  then  PV  =  0  (2) 

G  G 

if  -80  SAP  <0  then  PV  =  (  %-^  U  P  +  42.  8  (3) 

Li  Lj  oU.  O 

if  0  s  4P  then  PV  =  42.  8  mm  Hg  (4) 

G  G 

Note  that  42.8  mm  of  Hg  is  equivalent  to  2G  protection  since  the  cardio¬ 
vascular  system  gain  is  -?1.4.  The  simulation  model  is  presented  in 
block  diagram  form  in  Figure  4. 
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ACCELERATION  FORCE 


1 


Figure  3.  Standard  Pressure  Suit-Valve  Fill  Schedule 


Figure  4.  G-Suit  Model  for  Generation  of  Protection  Value  PV^ 
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2.3  THE  STRAINING  MODEL 

Human  tolerance  of  acceleration  may  be  improved  by  the  use  of  pro¬ 
tective  measures  and  devices.  Straining  and  grunting  in  tight  turns  was 
practiced  by  German  pilots  prior  to  World  War  II  to  improve  their 
tolerance.  This  technique  gradually  evolved  into  the  M-l  maneuver.  At 
the  same  time  it  was  discovered  that  bending  forward  also  improved  the 
tolerance  of  pilots  in  high  acceleration  maneuvers.  The  effects  ot  posture, 
restraint,  and  body  position  on  G  tolerance  have  been  well  worked  out  and 
have  been  reviewed  in  detail  (Davson,  1969;  Fraser,  1966;  and  Gaur  and 
Zuidema,  1961). 

Since  the  early  days  of  flying,  the  value  of  muscular  straining  as 
protection  against  the  effects  of  J-G  acceleration  has  been  recognized. 

Pilots  observed  that  blackout  or  grayout  could  be  postponed  if  they  grunted 
or  screamed  and  tensed  the  skeletal  musculature  during  acceleration.  The 
M-l  maneuver,  defined  as  muscular  straining  with  expiration  against  a 
partially  closed  glottis  (see  Appendix  C)  is  a  most  effective  voluntary  pro¬ 
tection  against  the  circulatory  effects  of  +G  acceleration  (Lindberg,  Sutterer 
Sutterer,  Marshall,  Headley,  and  Wood,  I960;  White,  1958).  The  physi¬ 
ological  basis  for  this  protection  has  been  ascribed  to  its  effect  on 
increasing  the  arterial  pressure  at  eye  level  during  acceleration.  As 
forced  expiration  is  instituted,  the  resulting  increase  in  intrathoracic 
pressure  is  transmitted  directly  to  the  aorta.  There  is  an  immediate 
increase  in  peripheral  arterial  pressure  which  is  only  slightly  less  than 
the  increase  in  intrathoracic  pressure.  This  maneuver  affords  up  to  2.4  G 
protection  against  blackout  (Burton,  et  al.,  1974).  The  circulatory 
mechanisms  described  were  identified  by  Rushmer  in  1946.  The  "pressure 
raising"  approach  has  theoretical  difficulties  in  that  blood  pressure  is 
controlled  within  rather  narrow  limits  by  a  number  of  effective  servo¬ 
mechanisms,  the  baroreceptor  system  being  the  best  known.  As  a  matter 
of  course  the  elevated  intrathoracic  pressure  can  result  in  a  decreased 
venous  return,  lowered  cardiac  output  and  therefore  a  lower  blood 
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pressure  at  the  eye.  The  carotid  sinus  reacts  to  the  initially  increased 
blood  pressure  and  also  attempts  regulation.  Both  of  these  effects  are 
minimized  by  a  cyclic  application  of  the  M- 1  or  L- 1  procedure  with  a 
repetition  period  of  4-5  seconds. 

The  skeletal  muscle  is  controlled  by  signals  which  are  transmitted 
to  selected  motor  units  of  a  muscle  through  the  motor  neurons.  The  force 
generated  by  the  muscle  is  the  result  of  both  the  frequency  of  firing  of 
motor  units  and  the  number  of  motor  units  which  are  recruited.  The 
muscle  tension  is  accompanied  by  electrical  signals  (EMG)  which  can  be 
detected  by  suitable  electrodes  on  the  skin  surface.  The  electrical  signal 
exhibits  the  characteristics  of  its  primary  source  in  that  it  represents  a 
weighted  sum  of  motor  unit  activations  (Shannon,  1975).  There  is  therefore 
a  correspondence  between  the  EMG  and  muscle  force  (Bigland  and 
Iippold,  1954). 

The  EMG  signal  can  be  considered  to  be  the  summation  of  a  large 
number  of  individual  potential  waveforms  each  resulting  from  the  firing  of 
a  motor  unit.  Thus  we  can  assume  that  the  EMG  signal  is  Gaussian  with 
a  mean  of  zero.  If  such  a  signal  is  the  input  to  a  full  wave  rectifier,  it  has 
been  shown  (Davenport  and  Root,  1958)  that  the  mean  output  is  proportional 
to  the  standard  deviation  of  the  input.  Furthermore,  the  standard  deviation 
has  been  shown  (Shannon,  1975)  to  be  the  appropriate  measure  of  the  EMG 
signal  to  use  for  control  purposes.  This  reasoning  leads  to  selection  of 
the  rectified,  smoothed  EMG  signal  for  determining  a  straining  protection 
value . 

The  simulator  model  requires  the  presence  of  two  signals  to  generate 
the  straining  protection  value,  PVS  .  The  two  signals  represent  the  prime 
factors  which  are  present  in  a  properly  executed  straining  maneuver.  The 
EMG  signal  is  generated  by  the  straining  subject  and  processed  to  provide 
an  intermediate  protection  value  PVm.  The  protection  afforded  by  muscle 
straining  is  then  modified  by  the  effects  ol  the  breathing  maneuver  ( M  -11. 
The  protection  factor  generated  by  the  breathing  pattern  is  critically 
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controlled  by  the  straining  and  breathing  rate.  An  improperly  performed 
M- 1  can  actually  reduce  protection  and  result  in  a  detrimental  pressure 
drop. 

We  assume  the  thoracic  cavity  to  be  represented  by  a  flexible 
cylinder  surrounded  by  a  muscle  girdle.  The  muscle  straining  signal 
representing  muscle  force  around  the  thoracic  volume  is  then  assumed  to 
bear  a  linear  relation  to  the  internal  pressure  increase  at  low  straining 
levels.  At  higher  levels  a  maximum  pressure  is  achieved  and  increases 
in  muscle  straining  are  no  longer  effective.  The  pressure  straining 
relationship  used  in  the  model  is  shown  in  Figure  5  as  a  curve  with  a 
saturation  level  related  to  the  maximum  pressure  rise  in  the  thoracic 
cavity. 

The  effectiveness  of  the  straining  maneuver  is  also  governed  by  the 

repetition  rate.  The  straining  protection  value  PVm  is  therefore  modified 

by  a  function  dependent  upon  maneuver  repetition  period  K  (t  ).  Experi- 

Z  iv 

mentally,  K  (t  )  is  maximum  at  3-5  seconds  according  to  Gillingham 


Figure  5.  Pressure-Straining  Relationship  Based 
on  EMG  Statistical  Properties 


(Appendix  C).  Shorter  time  periods  do  not  allow  sufficient  time  for  the 
internal  pressure  to  rise  and  longer  time  periods  involve  countering 
responses  due  to  baroreceptor  feedback  and  reduced  venous  return  to  the 
heart.  The  function  representing  the  variation  of  K  with  maneuver 
repetition  rates  is  shown  in  Figure  6.  The  straining  simulation  block 
diagram  is  shown  in  Figure  7. 


2.4  THE  VISUAL  FIELD  MODEL 

Rogers,  1978  discusses  the  structure  and  functioning  of  the  eye  as  it 
forms  the  basis  for  a  usable  model  of  the  visual  field  under  acceleration. 

The  effect  that  acceleration  has  on  the  visual  apparatus  is  observed  in 
terms  of  tunnel  vision,  greyout  and  blackout.  During  the  periods  of 
grayout  there  are  also  decreases  in  visual  acuity  and  brightness  contrast 
detection  ability.  Tunnel  vision  means  that  loss  of  visual  perception 
occurs  first  in  the  peripheral  areas  of  the  visual  field  and  progresses 
inward.  Although  there  are  many  factors  related  to  the  anatomy,  psychology 
and  physiology  of  the  human  which  are  responsible  for  these  effects,  the 
structure  of  the  eye  is  a  primary  factor,  and  it  provides  a  suitable  model 
for  the  effects  of  acceleration  on  the  visual  field. 

The  primary  cause  of  the  effects  of  accelerationoon  vision  is  believed 
to  be  the  reduction  of  blood  supply  to  the  rods,  cones  and  cells  in  the  retinal 
layers.  This  has  been  confirmed  by  various  inve stigators  (Lewis  and 
Duane,  1956;  Newsom,  Leverett  and  Kirkland,  1968).  The  eye  is  supplied 
by  the  opthalmic  artery  which  has  two  branches  in  the  eye  itself.  One 
branch,  the  retinal  artery,  enters  the  eye  through  the  optic  nerve  and  pro¬ 
vides  the  primary  source  of  blood  to  the  retinal  layers.  The  other  branch, 
the  ciliary  artery,  supplies  the  outer  layers  of  the  eye  and  the  ciliary 
processes.  Figure  8  shows  the  retinal  artery  entering  through  the  optic 
disc  and  branching  out  to  supply  the  entire  retina. 

The  structure  at  the  retina  is  idealized  as  a  network  of  cylindrical 
hydraulic  pipes  with  a  central  supply  at  the  entrance  of  the  optic  nerve, 
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Figure  6.  Effect  of  Repetition  Period  on  Straining  Protection  Value 


□ 


Figure  8.  Topography  of  the  Eyeball 

from  W.  J.  White  in  Reference  SB 


Figure  9.  If  a  steady  laminar  flow  is  assumed  through  a  cylindrical  tube 
the  flow,  Q,  will  vary  with  the  4th  power  of  the  radius  of  the  tube  and  the 
pressure  differential  AP.  Also  the  flow  will  vary  inversely  as  the  length  A 
of  the  tube  and  the  kinematic  viscosity  of  the  fluid  r,.  This  is  stated  in 
Poiseuille's  law  as: 
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The  blood  pressure  drop  across  the  cardiovascular  system  is  mostly 
across  the  arterioles  (Berne  and  Levy,  1972).  At  higher  levels  of  pressure 
the  retinal  arteries  supply  the  arterioles  with  sufficient  pressure  that  the 
entire  retina  is  perfused  with  a  continuous  flow  of  freshly  oxygenated  blood. 
As  the  pressure  at  the  retinal  artery  drops,  the  flow  must  decrease  in 
accordance  with  Poiseuille's  Law,  and  less  blood,  therefore  less  oxygen, 
is  available  to  replenish  the  retina. 

The  effective  pressure  (P  )  in  the  retinal  artery  is  the  difference 

K. 

between  supply  pressure  P  and  the  interocular  pressure  P  (a  pressure 

a.e  IO 

differential  of  10-20  mm  Hg  is  required  to  maintain  eyeball  rigidity). 


As  an  initial  approximation,  the  retinal  supply  network  is  assumed 
to  be  a  network  with  the  supply  pressure  distributed  linearly  across  the 
network.  That  is  at  any  point  along  the  supply. 

dP 

r 

— : —  =  constant 
dx 


where  x  is  the  distance  in  degrees  along  the  arteriolar  bed,  P^  is  the 
corresponding  differential  pressure  along  x,and  the  slope  is  a  function  of 
the  supply  pressure  at  the  optic  disc. 

The  family  of  curves  describing  the  pressure  gradient  across  the 
retina  is  a  series  of  straight  lines  with  slope 
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where  P_  represents  the  current  retinal  arterial  pressure  at  the  disc  and 
K 

K  is  the  outer  circulation  limit  of  the  retina  corresponding  with  the  ora 

T: 

seratta.  For  any  value  of  P  the  retinal  pressure  as  a  function  of  degrees 

2L  6 

from  the  disc  is  given  by; 
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In  a  small  blood  vessel  a  pressure  differential  of  about  15-20  mm  Hg 

across  the  vessel  wall  is  required  to  maintain  an  open  flow  channel  (Berne 

and  Levy,  1972),  When  this  transmural  pressure  is  not  maintained,  the  vessel 

collapses,  and  the  blood  supply  is  cut  off  to  all  points  beyond.  This  critical 

pressure  differential  is  now  labeled  P_  . 

CRIT 


The  threshold  value  for  any  x  (called  x  )  in  Figure  10  occurs  when 

±J 

P  falls  below  P  .  Solving  for  x  in  terms  of  P  at  the  P  : 

t*  KJ.  I  J_j  ae  CKJL 1 
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For  a  nominal  value  of  =  110°  representing  the  outer  circulation  limit, 
xL  is  given  by:  (let  PJO  =  20  mm  Hg,  PCRIT  =  !7  mm  Hg) 


110(P  -  37) 

_ ae 

P  -  20 
ae 


(9) 


The  quantity  x^  represents  arc  length  measured  from  the  center  of 
the  optic  disk.  The  optic  disk  is  located  on  a  visual  field  map  approximately 
10°  to  the  right  of  the  optical  axis  and  10°  below  the  horizontal  right  eye 
and  to  the  left  and  below  for  the  left  eye.  The  peripheral  edge  of  the  model 
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Figure  10,  Poiseuille's  Law  Applied  to  the  Retina 

(where  dP/dx  =  K.  for  Various  P  ) 
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is  determined  by  a  polar  projection  of  the  retinal  field.  The  coordinate 
origin  of  the  polar  map  is  taken  as  the  visual  axis  and  the  optic  disc  is 
taken  as  the  flow  center.  Defining  as  degrees  of  angular  location 
measured  counterclockwise  about  the  optical  axis  (Figure  11),  the  field 
contribution  of  the  right  eye  under  acceleration  is  then  given  by  the  solution 
of: 
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2r  ,r  cos  (d  +  ty)  +  r  , 
d  a 


(10) 


where  r  (degrees)  is  the  sensed  visual  limit  at  ,  r^  is  offset  distance 
(degrees)  of  the  disc  from  the  visual  axis  and  d  is  the  angular  offset 
(degrees)  of  the  disc  from  the  horizontal  axis.  A  similar  expression  is 
found  for  the  left  eye. 

The  visual  field  model  can  now  be  extended  to  a  binocular  field. 

Figure  12  is  a  binocular  field  map  showing  the  coincident  foveal  areas  as 
the  visual  center,  and  the  optic  disc  for  each  eye  in  a  manner  that  depicts 
the  observed  field  from  inside  out.  The  outer  lines  depict  the  outer  edges 
of  the  peripheral  vision.  Lowered  blood  pressure  supply  in  each  of  the 
modeled  retinas  causes  the  field  to  collapse  as  concentric  circles  with 
centers  at  each  of  the  optic  discs  as  shown  by  the  dashed  lines.  Thus  the 
visual  field  collapses  toward  a  somewhat  oval  shape  with  the  vertical  field 
having  a  smaller  visual  angle  than  the  horizontal  field.  This  shape  is  in 
agreement  with  Gillingham's  experimental  results  (Gillingham  and 
McNaughton,  1972). 

At  some  point  after  the  visual  field  has  collapsed, the  central  vision 
is  lost.  Gillingham,  et  al.,  1972,  shows  that  central  vision  remains  after 
the  peripheral  field  has  collapsed.  The  operating  model  is  therefore 
modified  to  allow  for  central  vision  by  incorporating  a  central  visual  cone 
which  is  responsive  to  the  eye  elevel  protection  value.  Mechanization  of 
the  model  utilizes  a  binary  decision  to  determine  presence  of  the  central 
field.  The  decision  point  is  chosen  at  a  nominal  10  mm  Hg  below  that  of 
the  peripheral  field  critical  point.  j 

] 

i 


22 


In  the  actual  mechanization  of  the  visual  field  model,  the  following 
data  are  used: 

visual  cone  for  foveal  vision  =  5°  (Gillingham,  et  al.,  1972) 

rd  =  10°  (for  both  left  anc*  ri8ht  eyes) 

=  .  175  radians 

d  =  10°  (right  eye)  measured  clockwise  from  the  right  horizontal 
visual  axis 
=  .  175  radians 

d  =  170°  (left  eye )  measured  as  for  right  eye 
=  2.  967  radians 


110  (P  -  37) 

ae  , 

*L  =  — p - 7^0 -  degrees  (Equation  9) 

ae 

=  1.920  (P  -  37)/P  -  20)  radians 

ae  ae 

P  in  this  equation  becomes  PV,  the  generated  protection 

3-6 

value,  which  is  equivalent  to  eye -level  blood  pressure  under 
acceleration. 


The  above  information  is  used  in  Equation  10  which  is  repeated 

below: 
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2r  r  cos  (d  +  )  -t- 
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(10) 


All  angles  in  degrees  are  converted  to  angles  in  radians  and  the  equation 
is  solved  for  "r",  which  is  the  visual  limit  at  angle  t  about  the  visual  axis. 

In  either  field: 


✓  2  2  2  2 

r  cos  (d  4  ijr)  -  (r  -  xT  )  (11) 

V4.  Cl  a  Li 

In  the  right  field: 

r^_  -  .  175  cos  (.  175  +  f) 

_  -  PV  -  3  7  2 

.031  cos  (.175  +  t)  -  .  031  +  3.686  (  )  (12) 
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In  the  left  field: 


r  =  .  175  cos  (2.967  +  ) 

l 

.031  cos  (2.967  +  *)  -  .  031  +  3.686  (  g— -  1  (13) 

Thus,  using  Equations  12  and  13,  the  sensed  visual  limit  may  be 
calculated  for  any  angle  \|r  and  for  any  protection  value  PV.  These  equations 
are  applied  only  to  the  point  where  all  peripheral  vision  has  collapsed  and 
only  the  central  (foveal)  vision  remains.  The  visual  cone  for  foveal  vision 
is  selected  as  5°  (Gillingham,  et  al.,  1972).  In  the  software,  the  collapse 
point  is  reached  at  a  value  of  PV  which  leaves  the  quantity  under  the  radical 
in  Equation  12  positive: 


-.031  +  3.686  ( 


PV -3  7  2 
PV -20 ' 


=  0 


which  gives. 


(14) 


PV  =  38.  72  (15) 

The  relationship  between  retinal  blood  perfusion,  and  the  ability  for 
contrast  detection  in  a  visual  scene  is  apparently  an  oxygen  transport 
phenomenon.  The  observed  visual  dimming  and  field  collapse  model  could 
possibly  be  extended  to  include  these  effects.  Chambers  and  Hitchcock,  1963, 
have  documented  the  relationship  between  contrast  detection  threshold  and 
the  G  field,  and  White  (1958),  has  documented  the  effect  of  lowered  visual 
acuity  due  to  reduced  oxygen  content  of  blood  in  the  retinal  artery.  Using 
the  CV  model  for  the  response  of  eye  level  blood  pressure  to  acceleration, 
the  contrast  and  acuity  results  could  be  added  to  the  visual  dimming  and 
field  collapse  model. 

In  the  present  simplified  mechanization,  it  has  been  assumed  that  the 
brightness  and  contrast,  at  any  given  point  in  the  visual  field,  are  linearly 
related  to  PV,  the  eye  level  blood  pressure  as  influenced  by  acceleration. 
When  PV  is  known,  the  Visual  Response  Algorithm  produces,  in  addition  to  the 
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visual  field  limit,  signals  appropriately  scaled  to  be  used  for  controlling 
both  brightness  and  contrast  at  any  point  in  the  visual  field. 

The  algorithm  presently  includes  control  of  a  center  and  two  side 
visual  TV  screens  with  separate  brightness  and  contrast  signals.  Collapse 
of  the  peripheral  vision  will  result  in  the  two  side  screens  being  dark 
while  collapse  of  the  foveal  vision  will  result  in  all  three  screens  being 
dark.  Prior  to  collapse,  the  screens  would  exhibit  a  linear  dimming. 

2.5  THE  E MG  SUBSYSTEM 

Figure  1  shows  the  EMG  subsystem  between  the  pilot  and  the  strain¬ 
ing  model.  The  task  of  the  EMG  subsystem  is  to  sense  the  raw  EMG  signal 
from  the  pilot's  skin,  amplify  it,  process  it;  and  present  it  to  the  digital 
computer. 

The  expanded  block  diagram  of  the  EMG  subsystem  is  shown  in 
Figure  13.  Three  electrodes  (one  is  a  reference  ground)  pick  up  the 
pilot's  EMG  signals.  The  bio-amplifiers  add  the  two  electrode  voltages 
and  provide  a  gain  of  100.  A  hi-pass  filter  eliminates  low  frequency  drift. 
Amplification  up  to  1000  then  readies  the  signal  for  full-wave  rectification. 
Finally,  a  two- stage  low- pass  filter  provides  a  smooth  usable  signal  for 
use  in  the  simulator  computer. 


AMP.  HI  -  PASS;T|  =  .15  L0-  R*SS;T2  =  .8 

Figure  13.  EMG  Subsystem  Block  Diagram 

Except  for  the  electrodes  and  the  bio-amplifiers,  the  EMG  sub¬ 
system,  as  shown  in  Figure  13,  was  patched  on  an  EAI  TR-20  Analog 
Computer.  The  patchboard  wiring  diagram  is  shown  in  Figure  14. 
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Figure  14.  TR-20  Wiring  Diagram  -  EMG  Subsystem 


2.6  SYSTEM  INTEGRATION 


The  required  dynamic  model  of  pilot  visual  response  to  G  can  now  be 
assembled  as  a  summation  of  the  separate  models.  The  completed  system 
algorithms  provide  the  means  for  implementation  in  aircraft  training  simula¬ 
tors.  The  integrated  system  is  shown  in  Figure  15. 


Pilot  commands  enter  the  simulator  aircraft  dynamics  ^2)  causing 

changes  in  the  aircraft  velocity  vector.  The  resulting  g-loading  G(t),  is  used 
to  drive  the  cardiovascular  model  .  In  addition  the  G(t)  signal  causes  the 
pilots  G-suit  to  inflate  according  to  a  preselected  schedule.  The  effect 

of  the  proper  suit  inflation  is  to  increase  the  value  of  the  protection  variable 
within  certain  limits.  The  sum  of  these  signals  is  added  with  a  computer 

generated  random  signal  related  to  individual  variations  and  with  a  set  value 
representing  the  eye-level,  nonaccelerating,  blood  pressure  (120  mm  Hg). 

The  resulting  PV  signal  drives  the  visual  field  model.  The  instrument  panel 
and  the  view  screen  are  then  dimmed  according  to  the  output  of  the  model  ©• 
The  pilot  senses  the  size  and  brightness  changes  and  has  the  option  of  reducing 
his  aircraft  maneuver  intensity  (g-loading)  or  increasing  his  g-tolerance  level 
by  performing  a  straining  maneuver  (M-l  or  L-l).  If  he  initiates  such  a  man¬ 
euver,  the  EMG  subsystem  processes  the  myoelectric  activity  from  his  skin 
surface,  providing  a  smoothed  EMG  signal  0  to  the  straining  model.  The 
straining  model  then  generates  a  protection  value,  PV  ,  ©  ,  which  is  depen¬ 
dent  on  the  straining  level  and  on  the  straining  interval.  This  PV^  is  added 
with  the  other  effects  and  thus  serves  to  increase  PV  and  to  increase  the 
visual  field  limits  and  brightness. 


C ircled  number s  refer  to  Figure  15. 


SECTION  3 


ALGORITHM  DEVELOPMENT  AND  DOCUMENTATION 

3.  1  INTRODUCTION 

Two  FORTRAN  algorithms  have  been  developed  from  the  physiolog¬ 
ical  models  presented  in  Section  2.  These  algorithms  are: 

1.  The  Dynamic  Visual  Response  Algorithm 

2.  The  EMG  Algorithm 

The  Dynamic  Visual  Response  Algorithm  consists  of  a  FORTRAN 
coded  implementation  of  the  Cardiovascular  Response  Model,  the  Protective 
Garment  Model,  and  the  Visual  Response  Model.  The  EMG  Algorithm 
consists  of  a  FORTRAN  coded  implementation  of  the  EMG  Straining  Model. 

A  description  of  each  algorithm  is  presented  along  with  its  flow  chart. 

In  addition,  a  program  listing  combining  the  two  algor ithm^,  and  the  program 
results  with  several  input  profiles^  are  presented. 

3.2  DYNAMIC  VISUAL  RESPONSE  ALGORITHM 

The  flow  chart  for  the  Dynamic  Visual  Response  Algorithm  appears 
in  Figure  16.  The  first  input  is  simply  the  acceleration  in  units  of  G's 
(G  -  9.81  m/sec  [Rogers,  1978]).  This  input  must  be  generated  from  the 
aircraft  dynamics.  This  algorithm  assumes  that  the  acceleration  is  readily 
available.  The  acceleration  is  used  to  determine  the  cardiovascular 
protection  variable  contribution,  PVC.  The  cardiovascular  response 
model  is  represented  by  a  linear  transfer  function  which  corresponds  to 
the  response  of  the  pilot's  blood  pressure  to  acceleration.  The  calculations 
used  in  the  transfer  function  are  implemented  with  a  recursion  relation. 

As  with  all  the  various  protection  variable  contributions  a  positive  increas¬ 
ing  value  indicates  an  increased  protection  or  tolerance  level. 
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Figure  16.  Dynamic  Visual  Response  Algorithm  (concluded) 
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The  second  input  is  the  difference  between  the  subject's  actual  G-suit 
pressure  and  the  "standard  curse  pressure"  (Rogers,  1978)  (as  calculated 
from  the  standard  pressure  suit  valve  fill  schedule).  This  difference, 
which  is  labeled  DELP,  is  used  to  calculate  another  protection  variable 
contribution,  PVG.  The  final  input  is  from  the  EMG  algorithm  and  consists 
of  a  single  protection  variable  contribution,  PVS.  This  input  is  explained 
more  completely  in  the  EMG  algorithm  section. 

After  all  the  protection  variable  contributions  have  been  calculated, 
they  are  summed  to  form  the  protection  variable,  PV .  This  protection 
variable  is  then  used  to  calculate  the  blood  pressure  at  eye  level,  PAE. 

The  variable  PAE  is  then  used  to  determine  the  size  of  the  visual  field. 

The  size  of  the  visual  field  is  then  used  to  determine  equivalent 
light  levels  for  center  and  side  CRT  screens.  Brightness  and  contrast 
signals  for  the  center  and  side  screens  are  then  calculated  by  scaling  these 
light  levels. 

This  algorithm  assumes  that  all  three  inputs  are  readily  accessible, 
and  the  protection  variables  are  positive.  This  algorithm  was  written  in 
Fortran  IV  and  was  tested  on  both  the  University  of  Dayton  Sperry-Univac 
70/7  Computer  and  the  Air  Force  Human  Resources  Laboratory  (AFHRL) 
Simulation  Training  and  Research  System  (STARS)  Xerox  Sigma  5  computer 
facility  at  Wright-Patterson  Air  Force  Base. 

3.3  EMG  ALGORITHM 

This  algorithm  is  designed  to  take  the  rectified  EMG  signal  from 
the  subject's  straining  muscles  and  generate  the  protection  variable  con¬ 
tribution  due  to  straining,  PVS.  The  variable  PVS  in  turn  will  be  used  in 
the  main  program  as  one  of  the  protection  variable  contributions.  In  the 
process  of  calculating  PVS,  this  program  also  calculates  the  EMG  variable 
PVM^  and  the  straining  effectiveness  value,  K2.  Both  these  variables, 

PVM  and  K2,  are  important  because  they  are  separate  functions  of  time. 


Their  product  determines  the  straining  protection  variable  contribution, 
PVS.  The  algorithm  flow  chart  appears  in  Figure  17. 

The  single  input  of  this  algorithm  is  the  rectified  Electromyographic 
(EMG)  signal  from  the  subject's  body.  The  EMC  signal  is  used  only  in 
the  determination  of  the  EMG  protection  variable  PVM.  Upon  the  initiation 
of  a  desired  level  of  straining  the  counter  variable  TIME  is  started.  After 
TIME  is  initiated  the  straining  effectiveness  value  K2  is  calculated  accordin 
to  the  current  value  of  TIME. 

The  single  output  of  this  algorithm  is  the  straining  protection 
variable  contribution,  PVS.  The  variable  PVS  is  calculated  with  the  above 
mentioned  variables  K2  and  PVM.  Also  a  scaling  variable,  R,  is  used  to 
keep  the  value  of  PVS  between  acceptable  levels.  Appropriate  limiters 
have  been  written  into  the  algorithm  to  keep  the  variable  PVS  positive. 

This  algorithm  was  written  in  Fortran  IV  and  was  tested  on  both  the 
University  of  Dayton  Sper r y- Un ivac  70/7  Computer  and  the  AFHRL  STARS 
facility. 

3.4  PROGRAM  CHECK-OUT 

These  two  algorithms  have  been  combined  together  into  a  single 
program  for  use  as  software  for  the  BIOCONAID  system.  The  listing  of 
this  program,  called  BIO,  is  shown  in  Appendix  B. 

Six  different  computer  generated  input  profiles  were  used  for 
program  check-out.  Three  of  these  had  combined  acceleration  and 
straining  inputs,  while  three  contained  acceleration  inputs  only.  The 
results  of  these  profiles,  in  terms  of  the  protection  variable,  PV 
(calculated  blood  pressure  at  eye  level),  are  presented  in  Figures  18  to 
23.  All  of  these  results  exhibit  a  delayed  loss  in  PV  due  to  acceleration 
and  a  delayed  increase  in  PV  due  to  straining.  In  addition,  these  results 
agree  with  the  models  presented  in  Section  2. 
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SECTION  4 

STRAINING  EXPERIMENTS  WITH  HUMAN  SUBJECTS 


4.  1  INTRODUCTION 

The  experiments  and  equipment  development  described  here  pertain 
to  the  EMG  subsystem  which  was  described  earlier  in  the  report  (Section 

2.  5).  These  experiments  were  done  to  accomplish  the  following  goals: 

1.  To  evaluate  possible  electrode  sites  on  the  pilot's  body  to  picK 
up  straining  signals. 

2.  To  show  that  a  reliable,  usable  signal  is  produced. 

3.  To  indicate  tolerance  of  the  EMG  subsystem  to  extraneous  body 
motions . 

4.  To  aid  in  the  development  of  electronic  filters  for  the  EMG 
subsystem. 

The  experiments  involved  recording  the  EMG  signals  produced  as  the 
subjects  strained  in  response  to  visual  commands.  The  seven  subjects 
were  faculty  members  and  students.  All  faculty  members  were  experi¬ 
enced  pilots.  The  work  was  accomplished  at  the  University  of  Dayton  and 
at  the  Air  Force  Human  Resources  Laboratory,  Wright-Patterson  Air 
Force  Base,  during  the  period  of  November  1978  to  February  1979. 

4.  2  EXPERIMENTAL  SET-UP  AND  PROCEDURE 

Figure  24  shows  the  experimental  setup.  The  equipment  used  and 
the  experimental  procedure  are  described  below. 

Equipment:  1.  Electrodes  and  wiring 

2.  Bio-signal  amplifiers  (Gain  =  100) 

3.  TR-20  Analog  (Hybrid)  Computer  (EAI)  to  provide: 

a.  Hi-Pass  Filter 

b.  Amplification 

c.  Signal  Rectifier  (full  wave) 

d.  Lo-Pass  Filter  for  smoothing 

e.  Voltmeter  -  used  for  visual  feedback  of  error  signal 
to  subject 

4.  Strip  Chart  Recorder 
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Figure  24.  Experimental  Setup 

(equipment  "on"  and  checked  out  earlieri 

1.  Prepare  subject.  Strap-on  electrodes.  Use  EMC  Lotion. 
Explain  M-l/L-l  straining  maneuvers.  Explain  that  his 
straining  will  make  error  needle  move. 

2.  Subject  practices  straining  and  observes  needle.  Operator 
inserts  practice  commands  and  subject  strains  to  bring 
needle  back  to  zero. 

3.  First  test  period  (Experiment  No.  1  ).  Electrodes  on 
back.  Subject  relaxed  and  motionless.  Command 
level  one  (3  volt)  is  inserted  and  subject  strains  to 
keep  error  zeroed.  After  about  five  seconds, 
command  is  removed  and  subject  relaxes.  Repeat, 
when  subject  ready,  using  command  level  two  ( (>  volt). 

4.  Secont  test  period  (Experiment  No.  2).  Repeat 
except  subject  periodically  raises  and  lowers  alter¬ 
nating  feet  (1-2  inches)  during  straining  periods. 

Change  electrodes  to  abdomen. 

5.  Third  test  period  (Experiment  No.  1).  Repeat  first. 

6.  Fourth  test  period  (Experiment  No.  4).  Repeat  second 
except  subject  raises  and  lower  alternating  arms 

(6-8  inches)  during  straining  period  instead  of 
raising  feet. 
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(5  minutes) 


Change  electrodes  to  buttocks  (seat). 


(2  minutes)  7.  Fifth  test  period  (Experiment  No.  5).  Repeat  first. 

(2  minutes)  8.  Sixth  test  period  (Experiment  No.  6).  Repeat  second 
except  extraneous  motion  is  turning  head  (about  45°) 
alternately  each  way  during  straining  period. 

(3  minutes)  9-  Remove  electrodes.  Equipment  off. 

(33  minutes  total  time) 

4.3  RESULTS 

Preliminary  results  are  presented  in  Figures  25-27  for  three  different 
subjects.  This  data  was  taken  with  an  early  technique  of  EMG  signal  pro¬ 
cessing.  The  upper  trace  is  the  commanded  signal,  which  the  subject 
attempts  to  match  with  his  smoothed  EMG  signal  (the  middle  trace).  The 
bottom  trace  is  the  raw  EMG  signal. 

Although  these  preliminary  results  gave  useful  raw  EMG  signals, 
the  smoothed  EMG  s  gnal  suffered  from  a  low  frequency  drift  and  excessive 
noise.  In  addition,  the  smoothed  EMG  signal  exhibited  an  excessive  time 
lag  compared  to  the  raw  EMG  signal. 

Several  changes  in  the  EMG  signal  processing  technique  overcame 
these  shortcomings  and  are  discussed  in  Section  4.4.  The  results  for 
experiments  using  five  subjects  conducted  at  the  University  of  Dayton  with 
the  final  signal  processing  technique  are  presented  in  Figures  28-  33. 

These  data  show  a  well-behaved  smoothed  EMG  signal,  free  oi  drift  or 
extraneous  noise. 

In  addition  a  number  of  runs  were  made  with  a  different  subject  in 
the  STARS  simulator  cockpit  at  AFHRL.  The  RIOCONAID  EMG  subsystem 
equipment  was  operated  at  the  STARS  facility  to  determine  whether  there 
would  be  an  excessive  amount  of  electronic  noise  interference  in  that 
environment.  These  results  are  presented  in  f  igures  34-40,  in  terms  of 
the  raw  EMG  signal  and  an  error  signal  (commanded  signal  minus 
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smoothed  EMG  output).  The  equipment  operated  satisfactorily  without 
shielded-cable s  and  amplifiers.  All  STARS  equipment,  including  the  TV 
view  screens,  was  operating  with  the  exception  of  the  cockpit  instruments 
(400  Hz). 

4.4  DISCUSSION 

Early  in  the  experiments  it  became  apparent  that  a  more  effective 
filter  was  needed  for  the  EMG  signal.  The  smoothed  EMG  signal  had  too 
much  noise  associated  with  it.  The  raw  EMG  signal  included  a  low 
frequency,  drift,  which  also  affected  the  smoothed  signal.  The  smoothed 
EMG  signal  also  exhibited  an  excessive  time  lag. 

Several  changes  were  made  to  the  EMG  signal  processing  technique 
to  improve  the  system.  A  high-pass  filter  was  added  to  remove  the  DC 
drift  from  the  incoming  EMG  signal.  The  low-pass  filter  (smoothing  net¬ 
work)  was  converted  to  a  two-stage  filter.  Time  constants  of  all  filters 
were  adjusted  for  near-optimum  response.  The  final  configuration  is  that 
shown  in  Figures  13  and  24.  Figures  28-40  show  a  smoothed  EMG  signal 
which  is  reliable,  free  from  drift  and  free  from  excessive  noise.  This 
smoothed  EMG  signal  shows  very  good  correlation  with  the  straining  (raw 
EMG)  signal  from  the  pilot. 

The  Experimental  Procedure  lists  Experiments  5  and  6  which  were 
to  be  done  with  electrodes  on  the  buttocks.  Only  Figures  25  and  26  show 
any  such  tests  being  accomplished.  The  reason  Experiments  5  and  6  were 
discontinued  is  that  consistently  much  better  results  were  obtained  with 
the  lattisimus  dorsi  (back)  muscles  and  with  the  external  oblique  (abdominal) 
muscles.  Figure  26  shows  a  comparison  between  Experiments  4  and  5. 

The  smoothed  EMG  response  to  the  3  volt  command  has  roughly  three 
times  the  magnitude  in  Experiment  4  (abdomen)  than  in  Experiment  5 
(buttocks).  Similarly,  the  raw  EMG  traces  showed  a  marked  difference. 
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#1  (Preliminary) 


Raw  EMG  Smoothed  EMG  Command  Raw  EMG  Smoothed  EMG  Command 


Another  problem  encountered  with  placing  electrodes  on  the  buttocks 
was  the  difficulty  in  strapping  them  on.  This  problem  would  lx*  eliminated 
with  the  adhesive  electrodes  used  in  later  tests. 

It  appears  the  great  advantage  of  seat  electrodes  .  mild  be  their 
possible  use  without  the  necessity  of  mounting  electrodes  on  the  pilot's 
skin.  Efforts  to  dev  elop  "seat  pan"  electrodes  should  l  e  encouraged.  If 
electrodes  must  be  mounted  on  the  pilot,  then  the  abdomen  or  back  vivos  a 
much  stronger  signal. 

Experiments  1  and  2  had  electrodes  on  the  back  while  Experiments  > 
and  4  had  electrodes  on  the  abdomen.  Examination  of  the  raw  EMC  traces 
in  Figures  28  through  33  indicates  either  site  is  satisfactory  but  the 
abdomen  usually  gives  a  slightly  stronger  sienal. 

Experiments  2  and  4  involved  having  the  subject  do  extraneous  bodv 
motions  while  straining.  Rhythmically  raising  and  lowering  feet  (Experiment 
2)  and  raising  and  lowering  arms  (Experiment  4)  were  the  specific  motions. 
Experiment  6  involved  rhythmically  turning  the  head  from  side  to  side. 

There  was  no  noticeable  effect  on  the  smoothed  EMG  signal  due  to  these 
motions.  However,  major  motions  (standing  up  or  bending  far  forward) 
would  saturate  the  signal  processing  equipment.  Also  touching  or  jostling 
the  electrodes  caused  a  large  response.  We  feel  that  the  normal  motions 
of  a  pilot  in  the  simulator  cockpit  will  not  adversely  affect  the  BIOCONAID 
system. 

Two  input  command  levels  (usually  3  and  6  volts)  were  used  to  show' 
that  the  pilot  can  control  his  level  of  straining  with  reasonable  accuracy. 
Figures  28  through  3?  show  the  smoothed  EMG  signal  approximating  the 
command  level. 

Adhesive  electrodes  were  obtained  in  time  to  be  used  in  part  of  the 
University  tests  on  one  subject  (Figure  33)  and  on  all  the  tests  in  the  STARS 
facility  at  AFHRL.  As  might  be  expected,  these  electrodes  maintained  a 
more  constant  skin  contact  then  the  strap-on  type  electrode.  Also,  they 
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Raw  EMG 


are  easier  to  mount,  may  be  used  several  times,  and  the  lotion  in  them 
may  be  replenished  if  dry.  However,  removal  of  these  adhesive  electrodes 
is  not  pleasant  since  they  pull  hairs  out. 

A  difficulty  with  amplifier  saturation  in  the  TR-20  was  initially 
encountered  at  the  STARS  facility.  However,  when  the  gain  was  reduced 
by  a  factor  of  10,  everything  operated  normally.  To  compensate  for  the 
reduced  gain,  command  signal  levels  were  also  reduced  to  0.3  and  0.6 
volts. 

Some  very  high  frequency  (2  MHz),  low  amplitude,  noise  pick-up  was 
observed  on  an  oscilloscope  connected  to  the  raw  EMG  signal.  This  may 
have  been  due  to  the  TV.  Since  it  did  not  interfere  with  the  BIOCONAID 
system  operation,  it  was  concluded  that  electronic  noise  pick-up  at  the 
STARS  facility  would  not  be  too  severe. 


SECTION  5 

CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  CONCLUSIONS 

The  BIOCONAID  System  offers  the  potential  to  provide  a  very 
realistic  simulation  of  visual  dimming  with  straining  during  large  accelera¬ 
tions.  The  Dynamic  Visual  Response  Algorithm  provides  a  simulated 
visual  response  in  agreement  with  experimental  observations  and  known 
physiology.  The  EMG  Algorithm  provides  an  input  to  the  visual  response 
which  is  also  in  agreement  with  experimental  observations  of  straining 
during  large  accelerations. 

The  EMG  subsystem  provides  a  smoothed,  reliable  straining  signal, 
not  adversely  affected  by  extraneous  limb  movements.  Either  abdomen 
(external  oblique)  or  the  back  (latissimus  dorsi)  muscles  are  suitable  as 
electrode  sites. 

All  parts  of  the  BIOCONAID  System  have  been  carefully  checked-out 
with  the  exception  of  the  Visual  Field  Model.  Equipment  limitations  at  the 
STARS  facility  prevented  the  full  implementation  of  the  Visual  Field  Model 
during  the  contract  period.  Thus  final  determination  of  system  effective¬ 
ness  must  await  complete  system  implementation  and  evaluation  in  an 
appropriate  simulator.  A  test  plan  for  such  implementation  and  evaluation 
is  contained  in  Appendix  A. 

5 .  2  RECOMMENDA  TION 

It  is  strongly  recommended  that  the  BIOCONAID  System  be  imple¬ 
mented,  demonstrated  and  evaluated  in  an  appropriate  Air  Force  simulator 
with  experienced  pilots.  Although  the  system  appears  to  provide  a  realistic 
simulation  of  visual  dimming  while  straining,  human  pilot  evaluations  are 
necessary  to  determine  system  effectiveness.  It  is  expected  that  such 


implementation  would  involve  fine  adjustment  of  some  algorithm  gains, 
thereby  enhancing  the  system. 


It  is  also  recommended  that  some  medical  research  be  conducted  into 
the  physiological  effects  of  sustained  straining  in  a  1-G  environment.  It 
would  be  extremely  useful  to  identify  the  level  of  straining  which  causes 
significant  detrimental  cardiovascular  effects,  as  well  as  the  level  of  strain 
ing  which  causes  negligible  detrimental  effects.  This  medical  research 
should  begin  with  a  literature  review  of  experimental  and  historical  data. 

If  this  review  does  not  yield  definitive  answers  on  the  above  straining  levels 
new  experiments  on  human  subjects  or  animals  should  be  conducted  to  yield 
the  desired  results.  This  information  is  necessary  for  the  safe  implementa 
tion  of  the  BIOCONAID  System  in  a  training  environment. 
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APPENDIX  A 


TEST  PLAN 
BIOC  ONAID 

(Bionic  Control  of  Acceleration  Induced  Dimming) 

SYSTEM  TEST  AND  EVALUATION 

A.  1  INTRODUCTION 

This  test  plan  pertains  to  the  performance  evaluation  of  the  final 
design  of  a  system  for  bionic  control  of  acceleration  induced  dimming 
(BIOCONAID).  Present  advanced  flight  simulators  provide  a  system  for 
dimming  the  pilot's  visual  scene  to  simulate  partial  blackout  when  the 
aircraft  normal  acceleration  exceeds  an  appropriate  level.  The  new 
system  goes  one  step  beyond  and  provides  for  control  of  the  dimming 
through  the  pilot's  muscular  straining  action  (M- 1  or  L- 1  maneuver). 

The  BIOCONAID  system  has  been  designed  by  the  University  of  Dayton 
under  Contract  F336  1 5-  77- C- 0080  with  the  Air  Force  Human  Resources 
Laboratory.  This  test  plan  has  been  prepared  to  fulfill  a  contractual 
requirement. 

The  plan  involves  installing  the  BIOCONAID  system  in  the  SAAC 
F-4  Simulator  at  Luke  AFB,  conducting  the  tests  using  experienced 
Air  Force  pilots,  removing  the  BIOCONAID  system  and  writing  a  report. 
The  intention  is  to  minimize  interference  with  normal  training  use  of  the 
simulator  during  the  evaluation.  A  seven  week  period  is  provided  for 
installation,  interfacing  and  check-out  followed  by  a  one-week  test 
period.  Disassembly  and  removal  will  be  accomplished  in  one  day. 

The  specific  time  schedule  for  the  tests  is  dependent  upon  the  date 
Air  Force  go-ahead  is  given  and  upon  other  Air  Force  instructions.  The 
University  of  Dayton  anticipates  this  evaluation  will  be  conducted  under 
an  additional  contractual  arrangement.  The  following  section  contains 

the  test  objectives. 
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A.  2  TEST  OBJECTIVES 


A.  2.  1  Objective  A 

To  evaluate  the  performance  of  the  final  design,  specifically, 
is  the  visual  scene  dimming  reliably  and  realistically  controlled  by  pilot 
straining?  To  answer  this  question,  the  system  must  be  observed  during 
long-term,  high-acceleration  turns.  In  addition,  however,  observations 
must  be  made  during  rapidly  changing  maneuvering  accelerations  as  well 
as  during  periods  of  low  dynamic  acceleration.  These  latter  observations 
will  be  for  the  purpose  of  determining  whether  the  final  design  is  free  from 
inappropriate  responses.  All  tests  will  be  done  with  and  without  the  electro¬ 
myographic  feedback  to  allow  comparisons  to  be  made. 

A.  2.  2  Objective  B 

To  conduct  other  tests  on  the  final  system.  These  "other 
tests"  include  further  evaluation  of  several  EMC  electrode  types  and  elec¬ 
trode  sites  and  further  evaluation  of  system  tolerance  to  extraneous  motion. 
These  tests  can  be  done  with  minimal  increase  in  test  time  beyond  the 
Objective  A  tests  and  they  should  be  very  significant  since  they  will  involve 
Air  Force  pilots  in  actual  simulated  flight.  Pilot  tolerance  for  the  electrodes 
is  the  most  important  question  to  be  answered.  Electrodes  to  be  included 
in  these  tests  are: 

•  strap-on,  reusable  many  times 

•  adhesive,  disposable 

•  electrode  belt 

•  seat  pan  electrodes  (if  available) 

Electrode  sites  included  are: 

•  abdomen  (external  oblique) 

•  back  (latissimus  dorsi) 

•  buttocks  (with  seat  pan  electrodes) 


A.  3  TEST  PROCEDURE 


A.  3.  1  General 

During  a  one-week  test  period  a  number  of  experienced  Air 
Force  pilots  will  act  as  test  subjects  to  perform  simulator  flights  aimed  at 
accomplishing  the  test  objectives.  The  planned  number  participating  is 
ten  (10),  although  more  would  be  welcome.  Each  pilot  will  be  briefed, 
electrodes  will  be  installed  and  a  short  practice  period  with  the  BIOCON- 
AID  system  will  be  provided.  Next,  the  actual  test  maneuvers  will  be 
accomplished.  These  maneuvers  consist  of  long-term,  high  accelera¬ 
tion  turns  in  both  directions,  short-term,  rapidly  changing  acceleration 
turns  or  aerobatics  and  normal  operations  near  one-G,  such  as  landing 
pattern.  Details  of  the  test  protocol  are  given  below.  Finally,  after  the 
flight,  the  electrodes  will  be  removed  and  the  pilot  will  be  debriefed. 

Each  pilot  will  provide  a  brief  written  observation  report  which  will  be 
included  in  the  final  test  report. 

The  actual  time -in -simulator  for  each  pilot  is  estimated  as 
30-45  minutes  and  the  total  time  per  pilot  is  estimated  as  1— 1-1/2  hours. 
The  time  of  day  for  the  tests  will  be  such  as  to  minimize  interfering  with 
the  normal  simulator  training  schedule.  Once  the  BIOCONAID  equipment 
is  installed,  it  can  be  used  or  not  used  as  desired.  Merely  placing  a  switch 
in  the  OFF  position  returns  the  simulator  to  its  original  configuration. 

A. 3. 2  Specific  Test  Protocol 

Time  (min)/Cum.  Activity 

13/13  Pilot  Briefing: 

-  Review  blackout  phenemena 

-  Review  M-l/L-1  maneuver 

-  Explain  BIOCONAID  system 

-  Explain  test  objectives 

-  Describe  test  procedure 

-  The  need  for  electrodes 


2/15 


Install  electrodes  according  to  the  following  table: 


10/25 


20/45 


No.  Pilots 

1 

1 

2 

3 

3 


lyp? 

Strap- on 
Strap- on 
Belt 

Paste- on 
Paste-on 


Location 

back 

abdomen 

waist 

back 

abdomen 


If  seat  pan  electrodes  are  available,  they  will  be 
used  with  all  pilots  (i.  e.  for  part  of  each  pilot-test 
period). 

Pilot  enters  and  starts  simulator,  climbs  to  5000 
ft  altitude  and  practices  turns  to  familiarize  him¬ 
self  with  BIOCONAID. 

Test  Period.  Pilot  does  the  following  maneuvers 
at  5000  ft,  Mach  .  8  on  instruction  from  controller 
(roll  level  after  each  turn): 

-  A  36 0°  turn  right  with  60°  bank  angle  (2G) 

-  Repeat  to  left  without  EMG  feedback 

-  A  360°  turn  right  with  70°  bank  angle  (2.  92  G) 
twist  body  and  head  to  look  over  right  shoulder 

-  Repeat  to  left  without  EMG  feedback,  twist 
body  and  head  to  look  over  left  shoulder 

-  A  360  turn  right  with  80  bank  angle  (5.  76G) 

-  Repeat  to  left  without  EMG  feedback 

-  A  series  of  alternating^right  and  left  turns  of 
45  with  bank  angle  70  ,  rolling  directly  from 
one  turn  into  the  next.  Eight  turns  will  be 
done  /last  four  without  EMG). 

-  A  Chandelle  right  followed  by  a  Chandelle 
left  (no  EMG) 

-  A  Lazy- Eight  (no  EMG  when  turning  left) 

-  A  rectangular  landing  pattern 


5/50 


Pilot  shuts  down  simulator,  deplanes,  and  electrodes 
are  removed. 


20/70 


Pilot  is  debriefed  and  gives  written  observations: 

-  Tolerance  for  electrodes. 

-  Was  energy  expenditure  to  avoid  blackout 
realistic  ? 

-  Compare  with  and  without  EMG  as  regards 
realism  of  simulated  flight. 

-  Was  the  onset  of  blackout  and  the  recovery 
through  straining  realistic? 

-  Was  the  beneficial  effect  of  the  G-suit  realis¬ 
tically  simulated? 

Tolerance  of  system  to  body  twisting. 

-  Further  comments  or  recommendations. 

A.  3.  3  Recordings 

Figure  41  shows  a  block  diagram  of  the  integrated  system 
including  pilot,  aircraft,  G-suit,  cardiovascular  model,  straining  model, 
and  visual  scene  model.  The  numbers  and  names  used  below  for  recorded 
variables  refer  to  Figure  41. 

Both  digital  and  analog  signals  will  be  recorded.  It  is 
anticipated  that  digital  variables  will  be  transferred  to  tape  and  printed  off¬ 
line.  The  analog  recorder  will  be  turned  on  and  off  by  one  of  the  operators 
as  the  test  progresses. 

Teletype  input  to  the  simulator  computer  will  result  in 
appropriate  printed  identification  of  data  pertaining  to  each  maneuver. 

When  a  new  maneuver  is  due,  the  teletype  input  is  changed  by  the  operator. 

Table  A-l  lists  variables  to  be  recorded  together  with  type 
(digital  or  analog)  and  frequency  of  observation. 

A.  4  TIME  SCHEDULE 

(Note:  T  =  First  test  day,  designated  by  AF) 

TIME  EVENT 

T-49  UD  personnel  (3)  visit  Luke  AFB  for  test  planning 

T-39  UD  ship  or  transport  equipment  to  Luke  AFB 
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TABLE  A-l 


TABLE  OF  RECORDINGS 

Each  Each 


Once 

Maneuver  Second  Continuous 

Date,  hour,  subject 

D 

X 

Electrode  Location 

D 

X 

Electrode  Type 

D 

X 

Maneuver 

D 

X 

Time 

D 

X 

Roll  Angle  (4> ) 

D 

X 

Heading  Angle  (t) 

D 

X 

G-level  (N^) 

D 

X 

Raw  Straining  Signal 

A 

X 

Processed  Strain  Signal 
(EMG) 

A 

X 

G-Suit  Protection  Value 

PVC 

D 

X 

Cardiovascular  Model 
Output,  PV^ 

D 

X 

Straining  Protection 
Value,  PVg 

D 

X 

Net  Protection  Value,  PV 

D 

X 

Dimming  Signal  (Output 

A 

X 

of  Visual  Scene  Model) 

Visual  Scene  Light 
Intensity 


a — " 


A 


X 


PROJECTION 


Figure  A- 1  System  Integration  Diagram 


T-37 


UD  engineer  /programmer  arrives  Luke  AFB,  begin  interfacing 
equipment  and  software 

T-27  Interfacing  complete 

T-  17  Piloted  system  checkout  starts  (one  AF  pilot  needed) 

T-  10  BIOCONAID  system  check-out  complete 

T- 5  Formal  demonstration  of  system  and  AF  approval  to  proceed 

with  test 

T  Testing  starts 

T+6  Testing  ends 

T+7  Equipment  removed,  simulator  restored  to  original  configuration 

T+67  Draft  copy  of  final  test  report  to  AFHRC 

A.  5  IMPACT  ON  SAAC  SIMULATOR  SCHEDULE 

It  is  expected  that  the  simulator  designated  for  these  tests  will  be 
already  fully  utilized  during  normal  duty  hours  for  training  activities. 
Consequently,  the  training  usage  should  retain  first  priority  and  the  test 
program  activities  would  proceed  on  a  "time  available"  basis.  Hope¬ 
fully,  evening  and  weekend  hours  would  be  available  regularly  during 
the  44  day  set-up  and  test  period.  It  is  estimated  that  an  average  of 
two  hours  operating  time  per  day  would  be  required. 

Care  will  be  taken  to  insure  that  installation  of  the  test  equipment 
and  software  is  done  in  such  a  way  that  it  does  not  affect  the  normal 
training  use  of  the  simulator.  This  means  that  simple  switches  and 
teletype  entries  will  disarm  all  test  provisions  and  leave  the  simulator 
in  its  original  operating  mode. 

Also,  the  test  provisions  will  be  installed  so  that  they  may  be 
rapidly  removed  leaving  the  simulator  in  its  original  configuration. 
Complete  removal  of  test  provisions  from  the  site  will  require  less 
than  one  day. 


81 


A.  6  AIR  FORCE  SUPPORT  REQUIREMENTS 


A.  6.1  Equipment 

SAAC  (F-4)  Simulator  {88  operating  hours) 

Includes:  -  digital  computer  capable  of  handling  software 

-  one  analog/digital  port 

-  light  dimming  capability 

Four  channel  strip  chart  recorder. 

A. 6. 2  Personnel 

1  Project  Officer  -  liaison  with  base  (part-time) 

1  Simulator  operator /maintenance  technician  -  to  operate 
simulator  (if  required),  monitor  test  installation  (part- 
time),  answer  questions 

10  (or  more)  voluntary  pilot  subjects  (1-1/2  hrs  each) 

1  pilot  for  "piloted  system  check-out"  prior  to  tests  (4  hrs) 

A.  6.  3  Work  Space 

Combination  office /workroom  near  simulator. 

Access  to  simulator, 

A. 6. 4  Security 

Secure  area  for  UD  equipment  (table  top  analog  computer, 
meters,  oscilloscope,  oscillator,  bio -amplifiers,  cables, 
tools). 

Security  passes  for  four  UD  project  workers. 

Car  pass  for  rental  car. 


A.  7  UNIVERSITY  OF  DAYTON  SUPPORT 

A.  7.1  Equipment  and  Supplies 

Oscillator 

Bio -Amplifiers  (6) 

Cables  and  wiring 
EMG  electrodes 
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A. 7. 2  Personnel 


2  engineers  (one  visit  per  week) 

2  technicians  or  students  (full-time  on  site) 

A.  7.  3  Transportation,  Housing  and  Food 

The  University  will  furnish  transportation  for  its  personnel 
and  equipment  as  well  as  housing  and  meals  for  its  personnel. 

A.  8  REPORT  REQUIREMENTS 

A  report  will  be  prepared  as  prescribed  by  Air  Force  instructions. 
A  draft  copy  will  be  furnished  to  AFHRL.  60  days  after  completion  of  the 
tests.  The  report  will  include  the  recorded  data,  the  data  analysis,  the 
pilot  comments,  conclusions  and  recommendations,  and  a  discussion. 
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APPENDIX  B 


00010  c 
00020 
00030 
00040 
00050 
0  0  06  0 
00070 
00080  C 
00090  C 
0 0100  C 

00110  c 
00120  C 
00130  C 
00140  C 
00150  C 
00160  C 
00170  C 
00180  C 
00190  C 
00200  C 
0021 0  C 
00220  C 
00230  C 
00240  C 
00250  C 
00260  C 
00270  C 
00280  C 
00290  C 
00300  C 
00310  C 
00320  C 
00330  C 
00340  C 
00350  C 
00360  C 
00370  C 
00380  C 
00390 
00400 


ALGORITHM  LISTINGS 


TFORSR.C .  E I  □  . DRV  ID  OUflM 

COMMON-'  I NPUTSvX  <6>  »  DELP-  PVS 
COMMON^OUTPUT^ERSn  BRC»  CTS» CTC 
DIMENSION  Y <3> 

RERL  N»NFM»K2,LLS!.LLC 
CB  =1.0 
CC  =  1 . 0 


INITIRLIZRTION 


G  =  GRIN  OF  CRRD 10 VASCULAR  TRANSFER  FUNCTION. 

R,B!.C>D-.E  =  TUSTIN  CO-EFFICIENTS  OF  TRANSFER  FUNCTION. 

F  =  OFFSET  DISTANCE  OF  THE  DPTIC  DISC  FROM  THE  CENTRAL 
VISUAL  AXIS  a  0  DEGREES  FOR  RIGHT  EYE) 

P  =  OFFSET  DISTANCE  OF  THE  OPTIC  DISC  FROM  THE  CENTRAL 
VISUAL  AXIS <110  DEGREES  FOR  THE  LEFT  EYE> 

SEAT  =  SEAT  ANGLE  CSTANDARD  IS  13  DEGREES> 

EYE  =  OFFSET  ANGLE  OF  EYE  FROM  PHYSIOLOGIC  Z-RXIS 
CTYPICALLY*  EYE  =  5  DEGREES) 

Z  =  CONVERSION  FACTOR  FOR  DEGREES  TO  RADIANS  CONVERSION. 
X <N>  =  G  =>  DIMENSIONLESS  MULTIPLE  OF  GRAVITY  CONSTANT. 

Y < 1 )  =  OUTPUT  OF  CARDIOVASCULAR  TRANSFER  FUNCTION. 

DELP  =  DIFFERENCE  BETWEEN  STANDARD  PRESSURE  SUIT  VALVE 
FILL  SCHEDULE  AND  THE  ACTUAL  SUIT  PRESSURE. 

PV  =  TOTAL  PROTECTION  VALUE. 

PVC  =  CARDIOVASCULAR  PROTECTION  VALUE. 

PVS  =  ELECTROMYOGRAPHIC  PROTECTION  VALUE. 

PVG  =  G-SUIT  PROTECTION  VALUE. 

PAE  =  CALCULATED  BLOOD  PRESSURE  AT  EYE  LEVEL. 

SI  =  POLAR  POSITION  ANGLE  IN  VISUALFIELD  MEASURED 

COUNTERCLOCKWISE  ABOUT  VISUAL  AXIS  FROM  RIGHT 
HORIZONTAL  AXIS.  <0. LE. SI . LT. 2»PI  RADIANS.) 

INITIALIZATION  FOR  DYNAMIC  VISUAL  RESPONSE  ALGORITHM 


G  =  -21.4 
A  =  -0.882609 


00410 
0  042  0 
00430 
00440 
00450 
00460 
00470 
00480 
00490 
00500 
00510 
00520 
00530 
00540 
00550 
00560 
00570  C 
00580  C 
00590  C 
00600  C 
00610  C 
00620 
00630 
00640 
00650 
00660  C 
00670  C 
00680  1 
00690  101 
00700 
00710  C 
00720  C 
00730 
00740  C 
00750  C 
00760  C 
00770  C 
00780  C 
00790  C 
00300  C 
00810  C 
00820 
00830  C 
00840  C 
00850  C 
00860 
00870 
00880 
00890  C 
00900  C 
00910  C 
00920  C 
00930  C 


B  =  1.875342 
C  =  -0.0947  *6 

D  =  0. 003634  ♦G 
E  =  0.09331  ♦S 
F  =  0.17453 
P  =  3.32 
SEAT  =  0.226893 
EYE  =  0. 0872665 
Z  =  57.2958 

x  a  >  =  o.o 

X<2>  =  0.0 
X  03}  =  0.0 
Y(l)  =  0.0 

Y  <2>  =  0.0 

Y  ('3>  =  0.0 
SI  =  0.0 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

INITIALIZATION  FOR  EMG  ALGORITHM 

N  =  0 
NNM  =  0 
NFM  =  5 
R  =  1.5 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

READ  <99* 1 01>  X  <1>  » EM6»  DELP 
FORMAT  <3  <F14.5> > 

xd)  =  xa>  -  i 

NNM  =  NNM  +  1 

EMG  ALGORITHM 

DETECT  STRAINING 

IF  f'EMG.  GT.  0.  7)  GO  TO  15 

NO  STRAINING  SO  LET  PVS  =0.0 

N  =  0 
PVS  =0.0 
GO  TO  21 

STRAINING  DETECTED  ==  START  COUNTER 
CALCULATE  AND  LIMIT  PVM 
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00940 

00950 

15 

00960 

C 

00970 

c 

00980 

c 

00990 
01000 
01010 
01  0£0 
01030 
01040 
01050 

c 

01060 

01070 

10 

01080 

01090 

11 

01100 

12 

OHIO 

C 

01120 

c 

01130 

01140 

01150 

20 

01160 

21 

01170 

C 

01180 

C 

01190 

C 

01200 

C 

01210 

C 

01220 

C 

01230 

C 

01240 

C 

01250 

C 

01260 

C 

01270 

C 

01280 

C 

01290 

C 

01300 

c: 

01310 

c 

01320 

c 

01330 

c 

01340 

01350 

01360 

c 

01370 

c 

01380 

c 

01390 

c 

01400 

c 

01410 

c 

01420 

c 

01430 

c 

01440 

c 

01450 

c 

PVH  *  <22. 39*EMG) -15.67 
IF  (PVM. GT. 51 . 5)  PVM  =  51.5 

DETERMINE  EFFECTIVENESS  DF  STRAINING  DUE  TO  TIME 
FACTOR  <K2) 

TIME  =N/NFM 

IF (TIME. LE. 2. 0>  GO  TO  10 

IF (TIME. LT. 4. 0>  GO  TO  11 

IF (TIME. LT. 5. 0>  GO  TO  12 

K2  =  <16.0  -  TIME) /I 1 . 0 

GO  TO  £0 

K2  =  0.5 

GO  TD  20 

K£=.25*TIME 

GO  TO  £0 

K2  =  1.0 

DETERMINE  AND  LIMIT  PVS 

PVS  =  K2  ♦  PVM  ♦  R 
IF  (PVS.LT.  0.  0>PVS  =  0.0 
N  =  N  +  1.0 
CONTINUE 


END  OF  EMG  ALGORITHM 


CARDIOVASCULAR  RESPONSE  ALGORITHM 
▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼▼ ▼▼▼▼ 

CALCULATE  PVC  FROM  CARDIOVASCULAR  TRANSFER  FUNCTION 


Y  (1)  =  A*Y<3)  +  B*Y<2)  +  C*X<3)  +  D^X (2)  +  E^X<1) 
PVC  =  1£0  +  Y<1) 


LIMIT  PVC  ==>  NOT  TO  BE  LESS  THAN  0 


IF  <  PVC.  LT.  0.0  )  PVC  =  0.0 
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01460  C 
01470  C 
01480  C 
01490  C 
01500  C 
01510  C 
01520  C 
01530  C 
01540  C 
01550 
01560 
01570 
01580  C 
01590  C 
01600  C 
01610  C 
01620  C 
01630  C 
01640  C 
01650 
01660  C 
01670  C 
01680  C 
01690  C 
01700  C 
01710  C 
01720  C 
01730 
01740  C 
01750  C 
01760  C 
01770 
01780 
01790 
01800 
01810 
01820  88 
01830 
01840  C 
01850  C 
01860  C 
01870  C 
01880  C 
01890  C 
01900  C 
01910  C 
01920  C 
01930  7 
01940 
01950 
01960 
01970 


G-SUIT  MODEL 


CftLCULftTE  PVG  ==>NOT  TO  EXCEED  42.8 


PVG=0. 

IF  <DELP.GT.-80>  PVG  =  <42.8*<DELP  +•  80>  /80>*X<1>  /7 
IF  CPVG. GT. 42. 9>  PVG  =  42.8 


CRLCULRTE  PV  —  TOTRL  PROTECTION  VRLUE 


PV  =  PVC  +  PVG  PVS 


CRLCULRTE  PRE 


PRE  *  PV  —  <0. 77*30. 48^C0SCSERT  -  EYE>*X<1>> 

IF  PRE  IS  LESS  THRN  37. 3 *  VISION  CDLLRPSES 

I F  (PV . GT . 37 . 3)  GO  TO  7 
LLS  =  0.0 
LLC  =  0.0 
WRITE <1 » 88) 

WRITE <98»87>  NNM» X <6> » EMG> PV 
FORMRT  CSX  *  'VISION  COLLRPSEO 
GO  TO  32 


CRLCULRTE  SENSED  VISURL  RNGLE  <RRDEG*RLDEG> 
WE  CRLCULRTE  THE  VISURL  LIMIT»RRDEG*RT  RNY 
RNGLE > SI  9  USING  ONLY  THE  RIGHT  EYE. 


W=1 . 91986*<PV-37> ✓ <PV-20> 

PI  =  3.14159 

IF  CSI.LT.PI/2.0>  SY  =  SI 
IF  CSI.GE.PI/2.0)  SY  =  PI-SI 
IF  <SI.GE.3.^PI/2.0>  SY  =  -2.«PI+SI 
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01980  C 
01990  C 
02000  C 
02010 
02020 
02030 
02040  C 
02050  C 
02060  C 
02070  C 
02080  C 
02090  C 
02100  C 
02110  C 
02120 
02130 
02140 
02150 
02160 
02170 
02180 
02190  32 
02200 
02210 
02220 
02230 
02240  C 
02250  C 
02260  C 
02270  € 
02280  C 
02290  C 
02300  C 
02310  8 
02320 
02330 
02340  87 
02350  33 
02360 
02370 
02380 
02390 
02400 
02410  C 
02420  C 
02430  C 
02440 
02450 
02460 


NOW  -PI'2. .LE.SY.LE.+PI/2. 

AA=F*COS<. 17453+SY> 

BB=SQRT<.  031*<CDS<.17453+SY>>**2-.  031+W**2> 

RIP.  =  AA  +  BB 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

CHOOSE  CORRECT  SENSED  VISUAL  ANGLE  AND  CONVERT 
TO  DEGREES 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

RRDEG  =  Z  ♦  R1R 

LLS  =  <PV-110.>/4.5+10. 

LLC  =  <CPV-85.  >  ••"'5.  0+1  0. 

IFCLLS.GE. 10. >  LLS  =  10. 

IF  CLLC.GE. 1 0. >  LLC  =  10. 

IF CLLS.LE. 0.>  LLS  =  0. 

IFCLLC.LE.  0.>  LLC=0. 

BRS=LLS»CB 

BRC  =  LLC  ♦  CB 

CTS  =  LLS  ♦  CC 

CTC  =  LLC  ♦  CC 

IF  <LLC. LE. 0. >  LLC  =  0. 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

_  OUTPUT  SENSED  VISUAL  ANGLE  AND  UPDATE  ARRAYS 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

X<6>  =  X<1)  +  1 

PRINT  X C6>  * EMG* PV» LLS> LLC 
WRITE  <98>87> NNM *  X  <6>  *  ENG » PV 
FORMAT <12* 3F9.3> 

X  <3>  =  X  <2> 

x  <2>  =  xa> 
xa>  =  o.o 

Y  <3>  =  Y  <2> 

Y  <£>  =  Y(l) 

Yd)  =  0.0 

RUN  FINISHED  ? 

I F  CNNM . NE . 99>  GO  TO  1 

STOP 

END 


APPENDIX  C 

THE  M-I  AND  L-l  STRAINING  MANEUVERS 

C.l  M-I  MANEUVER 

"Pilots  commonly  refer  to  the  M-l  Maneuver  as  the  'grunt'  maneuver 
since  it  approximates  the  physical  effort  required  to  lift  a  heavy  weight. 

The  M-l  maneuver  consists  of  pulling  the  head  down  between  the  shoulders, 
slowly  and  forcefully  exhaling  through  a  partially  closed  glottis,  and  simul¬ 
taneously  tensing  all  skeletal  muscles.  Pulling  the  head  downward  gives 
some  degree  of  postural  protection  (shortens  the  vertical  head-heart 
distance):  intrathoracic  pressure  is  increased  by  strong  muscular  expiratory 
efforts  against  a  partially  closed  glottis;  and  the  contraction  of  abdominal 
and  peripheral  muscles  raises  the  diaphragm  and  externally  compresses 
capacitance  vessels.  For  long-duration  G  exposures,  the  maneuver  must 
be  repeated  every  4  or  5  seconds.  When  properly  executed,  the  exhalation 
phase  of  the  M-l  results  in  an  intrathoracic  pressure  of  50  to  100  mm  Hg 
which  raises  the  arterial  blood  pressure  at  head  level  and  thereby  increases 
+G  tolerance  at  least  1.5  G.  The  inspiratory  phase  of  the  M-l  maneuver 
must  (be  done  as  rapidly  as  possible,  sic.  ).  '' 

C.2  L- 1  MANEUVER 

"The  L-l  maneuver  is  similar  to  the  M-l  maneuver  except  the  air¬ 
crew  member  forcefully  attempts  to  exhale  against  a  completely  closed 
glottis  while  tensing  all  skeletal  muscles.  Using  either  maneuver  the  pilot 
obtains  equal  protection,  i.e.  1.5  G  greater  than  relaxed  blackout  level  with 
or  without  the  anti-G  suit.  In  a  1972  study,  subjects  wearing  anti-G  suits 
and  performing  either  the  M-l  or  L-l  straining  maneuver  were  able  to 
maintain  adequate  vision  during  centrifuge  exposure  of  +9G  for  45  seconds. 
Higher  and  longer  runs  have  not  been  attempted.  However,  it  is  important 
to  note  (a  word  of  caution)  that  forcefully  exhaling  against  a  closed  glottis 
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without  vigorous  skeletal  muscular  tensing  (Valsalva  maneuver)  can  reduce 
+G  tolerance  and  lead  to  an  episode  of  unconsciousness  at  extremely  low 
levels.  Therefore,  instruction  and  training  on  the  proper  method  of 
performing  these  straining  maneuvers  is  essential." 

(Source:  Gillingham,  K.  K.  and  Krutz,  USAFSAM  AR  10-74,  Brooks  Air 
Force  Base,  Texas. ) 
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APPENDIX  D 


THE  SOLUTION  ALGORITHM  FOR  THE  CARDIOVASCULAR  MODEL 


A  digital  algorithm  is  used  to  represent  the  cardiovascular  transfer 
function. 


,D.l) 

1  +  3. 23s  +  5.  17s 

This  algorithm  employs  Tustin's  method^  for  computational  efficiency. 
For  a  given  recursion  period  T,  determined  by  the  simulation  computer 
frame  rate,  the  output  y(k)  is  given  by  the  following  equation. 


y(k)  =  -y(k-2) 


T2  -  6.46T  +  20.68 
T2  +  6.46T  +  20.  68 


y(k-l)  ~ 


2T 


41. 36 


T  +  6.46T  +  20.  68 


tx(k-2>V  -10-^T  -^(k-n 


2T 


T  +  6.46T  +  20.68 


T  +  6.46T  +  20.  68 


+  x(k) 


T  4  10.26T 
T**  +  6.46T  4  20.68 


(D.2) 


Here,  y(k)  represents  the  output  PVc  at  time  t^,  y(k-l)  is  PVc  at  tk_j  -  t^-T 
and  x(k)  is  G(t,  ).  A  complete  listing  of  coefficients  for  different  values  of 

K 

T  is  given  in  Table  D-l.  The  coefficients  are  derived  with  a  transfer 
function  gain  equal  to  1 . 

A  hypothetical  aircraft  maneuvering  profile  is  shown  in  Figure  D-l. 

The  predicted  eye  level  protection  value  (PV)  due  to  the  cardiovascular 
response  model  is  shown  in  Figure  D-2.  The  PV  values  are  derived  using 
T  =  .  02  and  the  hypothetical  profile  as  a  driving  function. 
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TABLE  D-l 

lasting  of  Coefficients  for  Pvc  Simulation 


T 

A 

B 

C 

D 

E 

0.  02 

0.9876 

-1.9875 

-0.0102 

0. 0000 

0. 0102 

0.04 

0.9753 

-1.9750 

-0. 0202 

0. 0002 

0. 0204 

0.  06 

0.9632 

-  1.9625 

-0.0301 

0. 0003 

0. 0304 

0.08 

0.9513 

-1.9501 

-0.0398 

0.0006 

0. 0404 

0.  10 

0.9394 

-1.9376 

-0.0493 

0. 0009 

0. 0502 

0.  12 

0. 9278 

-1.9251 

-0.0587 

0. 0013 

0.0600 

0.  14 

0.9163 

-1.9127 

-0. 0679 

0. 0018 

0. 0697 

0.  16 

0. 9049 

-1.9002 

-0. 0770 

0.0024 

0.  0793 

0.  18 

0.8937 

-1. 8878 

-0.0859 

0.0030 

0. 0889 

0.20 

0. 8826 

-1.8753 

-0.0947 

0. 0036 

0. 0983 

0.  22 

0.8717 

-1.8629 

-0. 1033 

0. 0044 

0.  1077 

0.  24 

0.8609 

-1. 8505 

-0. 1118 

0. 0052 

0. 1169 

0.26 

0.8502 

-1. 8382 

-0. 1201 

0. 0060 

0. 1261 

0.  28 

0. 8397 

-1. 8258 

-0.  1283 

0. 0069 

0. 1352 

0.  30 

0. 8293 

-1.8135 

-0. 1363 

0.0079 

0. 1443 

0.  32 

0. 8191 

-1. 8011 

-0.  1442 

0. 0090 

0. 1532 

0.  34 

0. 8089 

-1.7888 

-0.  1520 

0.0101 

0. 1621 

0.  36 

0. 799 0 

-1. 7766 

-0.  1597 

0. 01  12 

0. 1709 

0.  38 

0.7891 

-1.7643 

-0. 1672 

0.0124 

0. 1796 

0.40 

0.7794 

-1.7521 

-0. 1745 

0. 0137 

0. 1882 

0.42 

0.7698 

-1.7398 

-0.  1818 

0.0150 

0.  1967 

0.44 

0. 7603 

-1.7276 

-0. 1889 

0. 0163 

0. 2052 

0.46 

0. 7509 

-1. 7155 

-0. 1959 

0.  01 77 

0.  2136 

0.48 

0. 7417 

-1. 7033 

-0.  2027 

0. 0192 

0. 2219 

0.  50 

0. 7326 

-1.6912 

-0. 2094 

0. 0207 

0. 2301 

y(k)  =  -y(k  - 

■2)A  -  y(k-  1  )B 

+  xf k-21C  -* 

x(k-  1  ID  " 

x(klE 
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Figure  D“1  Hypothetical  Aircraft  Maneuvering  Profile 


